Polynomial Interpretation of Multipole Vectors 



Gabriel Katz 
Bennington College, Bennington VT 05201 USj^ 



o 
o 

(N 

C 
3 



Jeff Weeks 
15 Farmer Street, Canton NY 13617 USM 
(Dated: February 2, 2008) 

Copi, Huterer, Starkman and Schwarz introduced multipole vectors in a tensor context and used 
them to demonstrate that the first-year WMAP quadrupole and octopole planes align at roughly 
the 99.9% confidence level. In the present article the language of polynomials provides a new and 
independent derivation of the multipole vector concept. Bezout's Theorem supports an elementary 
proof that the multipole vectors exist and are unique (up to rescaling). The constructive nature 
of the proof leads to a fast, practical algorithm for computing multipole vectors. We illustrate the 
algorithm by finding exact solutions for some simple toy examples, and numerical solutions for the 
first-year WMAP quadrupole and octopole. We then apply our algorithm to Monte Carlo skies to 
independently re-confirm the estimate that the WMAP quadrupole and octopole planes align at the 
99.9% level. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 



The first-year WMAP data Q reveal a somewhat 
planar octopole which approximately aligns with the 
quadrupole |2|. More recent studies confirm these con- 
clusions at roughly the 99.9% level |9( while revealing 
mysterious alignments with the ecliptic plane 0, sug- 
gesting either a hitherto unknown solar system effect on 
the microwave background or an error in the collection 
and/or processing of the WMAP data. Other researchers 
find the I = 4 multipole is generic, the I = 5 multipole is 
unusually non-planar at the 99.8% level, and the I = 6 
multipole is unusually planar at the 98.6% level Q. No 
explanation is yet known for these strange results. 

Multipole vectors provide a convenient means to quan- 
tify the planarity of a given multipole, as well as to 
compare the alignment of two different multipoles 
The present authors, coming from a background in pure 
mathematics, were unable to decipher the formalism and 
terminology of Ref. p| and chose instead to re-create the 
multipole vector concept from scratch. The real-valued 
spherical harmonics of order I are precisely the homoge- 
neous harmonic polynomials of degree I in the variables 
x, y and z (for example Y 2 ° is the polynomial x 2 +y 2 — 2z 2 , 
up to normalization), so the present authors sought to 
understand the multipole vectors of Copi, Huterer and 
Starkman (CHS) from a polynomial point of view. 

Translated to the language of polynomials, CHS's mo- 
tivating goal (see Eqn. (10) of was to factor every 
homogeneous harmonic polynomial P of degree I into a 



product of linear factors 

P(x,y,z) = X ■ {a\x + byy + c\z) 

■ (a 2 x + b 2 y + c 2 z) 

■ (a e x + b e y + cez). 



(1) 



Such a factorization is of course impossible in general, 
as CHS implicitly acknowledge by their introduction of 
suitable error terms. In the language of polynomials the 
correct statement of the theorem is 

Theorem 1. Every homogeneous polynomial P of degree 
£ in x, y and z may be written as 

P(x,y,z) — A • (dxx + biy + c\z) 

■ (a 2 x + b 2 y + c 2 z) 

■ (a e x + b e y + cgz) 
+ {x 2 + y 2 + z 2 ) ■ R, 



(2) 
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where the remainder term R is a homogeneous polyno- 
mial of degree 1 — 2. The decomposition is unique up to 
reordering and rescaling the linear factors. 

Notes: (a) Theorem 1 lives entirely in the realm of real 
polynomials: the coefficients of P, R and all the linear 
factors aiX+biy+CiZ are assumed to be real, (b) Theorem 
1 does not require the polynomial P to be harmonic. 

In cosmological applications we are interested only in 
the value of the polynomial on the unit sphere S 2 ; we 
ignore its value on the rest of Euclidean 3-space. On the 
unit sphere the factor x 2 +y 2 +z 2 is identically 1, so in this 
case Theorem 1 says that any homogeneous polynomial 
P may be written as a product of linear factors \{a\x + 
b\y + c\z) ■ ■ ■ (agx + bey + cez) plus a remainder term R of 
lower degree. Applying this reasoning recursively gives 
the easy 
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Corollary 2. When restricted to the unit sphere, every 
polynomial P of degree £ in x, y and z may be written as 
P(x,y,z) = 

X e ■ (ag t ix + b^iy + c iA z) ■ (a^ 2 x + bt l2 y + C(.,%z) 
■ ■ ■ (a £ jx + bt,iy + ct,iz) 

+ ... 

+ A 2 • (a 2 ,ix + b 2 ^y + c 2 ,\z) ■ (a 2t2 x + 62,22/ + c 2,2z) 
+ Ai • (ax.jx + b 1A y + c lfl z) 

+ A . (3) 

The decomposition is unique up to reordering and reseat- 
ing the linear factors within each term. 

Note: Corollary 2 does not require the polynomial P to 
be either homogeneous or harmonic. 

Proof of Corollary 2. Write P as a sum of homogeneous 
terms P = Pi+Pi-i+- ■ -+Pi+P . First apply Theorem 1 
to the highest order term Pg, yielding a factorization Xi ■ 
(ae^x + bi^y + a^z) ■ ■ ■ (agjx + bejy + ci^z) along with a 
remainder term Rg- 2 of homogeneous degree £ — 2. (The 
factor x 2 + y 2 + z 2 may be ignored on the unit sphere.) 
Fold Ri- 2 in with Pi- 2 , and proceed recursively, applying 
Theorem 1 to P1-1, then Pi- 2 , and so on. 

To prove uniqueness, consider the even and odd parts 
of P separately. That is, write P = P C vcn + -fodd, where 
-Poven contains all the even-powered terms and P dd con- 
tains all the odd-powered terms. Say we have two poten- 
tially different decompositions for the even part 

^cvon = ir + ir_ 2 + n £ _4 + --- + n 

= IL' e +R' e _ 2 +W e _ i + --.+n' . (4) 

where each IT is the i th term in a decomposition and 
where the leading index will be £ or £ — 1 according to 
whether £ is even or odd. To make these decompositions 
homogeneous, multiply through by appropriate powers of 
Q = x 2 + y 2 + z 2 , 

P CV on = IT + QTle-2 + Q 2 n £ _4 + • • • + Q e/2 U 
= n' e + Qn' e _ 2 + Q 2 n' e _ 4 + ■■■ + Q e/2 n' . (5) 

This does not affect the value of P on the unit sphere, 
because Q — 1 there. The uniqueness part of Theorem 1 
implies that the leading order terms IT and W e must be 
equal. So subtract off those leading terms, divide through 
by Q, and apply Theorem 1 again to conclude Tle- 2 — 
H' e _ 2 - Continue recursively to finally reach IIo = IIq. 
The same argument then proves that the odd part of P 
has a unique decomposition as well. Q.E.D. 

II. PROOF OF THE MAIN THEOREM 

Even though the statement of Theorem 1 lives wholly 
in the world of real polynomials, its proof will dive deeply 
into the world of complex polynomials. So let the vari- 
ables x, y and z range over the complex numbers, while 



insisting that the coefficients of the polynomial P remain 
real. Because P has homogeneous degree £, whenever one 
point (xq, yo, zq) satisfies P(x,y,z) = 0, every nonzero 
constant multiple (axo, ayo, azo) satisfies it as well. Thus 
the equation P — is well defined on each equivalence 
class of points {a(a;o, yo, zq) \ a G C — {0}}. In other 
words, the complex curve P = is well defined on the 
complex projective plane CP 2 , which is the quotient of 
C 3 — {(0, 0, 0)} under the equivalence relation (xo, yo, Zq) 
~ ct(xo, yo, z o)- This leads us into the realm of algebraic 
geometry and puts its powerful tools at our disposal. 
The most useful tool for our purposes is 

Bezout's Theorem. If P and Q are homogeneous poly- 
nomials of degree m and n, respectively, then the curves 
P = and Q = intersect in CP 2 

• in exactly m ■ n points, counted with multiplicity, if 
P and Q share no common factor, or 

• in infinitely many points, if P and Q do share a 
common factor. 

For an elementary exposition of Bezout's Theorem, see 
□ 

In the present case, the only way the polynomial P 
may share a factor with the irreducible polynomial 
Q(x, y, z) = x 2 + y 2 + z 2 is for P to contain Q as a fac- 
tor, in which case Theorem 1 is trivially satisfied (take 
A = 0). So henceforth assume P does not contain Q as 
a factor. Bezout's Theorem now tells us that the degree 
£ complex curve P(x,y,z) = intersects the quadratic 
curve Q(x,y,z) — in exactly 2£ points, counted with 
multiplicities. None of the intersection points may be 
purely real, because real values cannot possibly satisfy 
x 2 + y 2 + z 2 = - recall that the definition of CP 2 ex- 
plicitly excludes (0,0,0). Furthermore, because P and 
Q both have real coefficients, whenever (xo, yo, zq) lies 
in the intersection P = Q = 0, its complex conjugate 
(jxo, yb~,~zb~) must lie there too. So the 2£ points of inter- 
section consist of £ pairs of non-real complex conjugates 

{Pi,¥T,---,Pi,Pl}- 

We claim that each pair {pi,pl} determines a unique 
line aiX + b.^y + CiZ = with real coefficients. The proof 
is easy. The conjugate pair {pi,jh} lies on the line aix + 
biJJ + ^z = if and only if the real and imaginary parts 
satisfy the following two totally real equations 

ai Rep^ x + bi Rcp hy + c t Rep ltZ = 

at Imp,j + hi lmp^ y + a lmp hZ = 0. (6) 

Geometrically those two equations represent planes in 
R 3 . If the coefficient vectors (Repi, x , Rep,;^, Kepi }Z ) and 
(Impi^ x , Impi^y, Impi >z ) are non-collinear, then the two 
planes are distinct and their intersection, which defines 
the solution set for (a^, bi,Ci), is a line through the origin 
in R 3 . In other words, the line aiX + biy + c^z = is 
unique. Normalize the coefficients to unit length, i.e. 
a 2 + b 2 + cf = 1, and the only remaining ambiguity is an 
overall factor of ±1. 
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But what if the coefficient vectors Repi = (Repi !X , 
Rep i>y , Rep hZ ) andlirift = (Im.p i<x , lm.p i>y , lmp hZ ) had 
been collinear? In this case the line aix + biy + CiZ = 
would be ill-defined. Fortunately this case does not 
arise. For if Ivapi were proportional to Repi, say Im pi = 
(3 Re pi, then the point pi, as an element of CP 2 , could 
be rewritten as a scalar multiple pi ~ \~^a Pi = R e Pii 
showing that pi is totally real. In other words, pi would 
lie in RP 2 C CP 2 . In particular, pi would be its own 
complex conjugate, and we can hardly expect a single 
point pi = pi to determine a unique line. Fortunately this 
case cannot occur, because Q(x, y,z) = x 2 + y 2 + z 2 = Q 
admits no real solutions. 

So let Li denote the unique line ciiX + biy + CiZ = 
containing the conjugate pair {pi,pt\. More precisely, 
let Li — ciiX + biy + CiZ be the unique (modulo rescaling) 
real linear polynomial whose roots include both pi and 
Pi. The desired decomposition J5J) becomes 

P = \L X L 2 ---L t + QR. (7) 

To prove that this equality holds, we again turn to 
Bezout's Theorem. First recall that the complex curve 
P = intersects the complex curve Q = in precisely 
the 2£ points {pi,Pi, ■ ■ ■ , PhPl}- By construction, the 
product curve L1L2 ■ ■ ■ Lg = also intersects Q = in 
those same 2£ points, and by Bezout's Theorem there are 
no other points of intersection. Now pick any other point 
q S {Q = 0} and define A to be the ratio 

A = — - — - — -—7— — -. (8) 

Li{q)L 2 {q) ■ ■ ■ Lg(q) 

Write a new polynomial 

F = P-\L 1 L 2 ---L l . (9) 

This new polynomial F has degree £, yet has zeros at the 
2£+l distinct points { q, p\,pi, . . . ,Pi,pg} C Q. In other 
words, the complex curve F = intersects the complex 
curve Q = at (at least) 2£ + 1 distinct points. By 
Bezout's Theorem the polynomials F and Q must share 
a common factor; because Q is irreducible the common 
factor must perforce be Q itself. Thus we may factor F 
as 

F = Q ■ R (10) 

for some remainder term R. Combining and 
yields the desired decomposition JJJ. 

Let us now prove that A is real. In light of the fac- 
torization (10), the polynomial F is clearly zero on the 
whole complex curve Q = 0. In particular, for the point 
q chosen earlier, 

F(q) = F(q) = 0. (11) 

On the one hand 

F(q) = P(q) - XL 1 (q)L 2 (q) ■ ■ ■ L e (q). (12) 



On the other hand, because P and the Li all have real 
coefficients, 

= P(q) - \L x {q)L 2 {q) ■ ■ ■ L e (q) (13) 

Comparing (|llfl . 1|12|) and i|13|) . and recalling that q was 
chosen to ensure Li(q) ^ 0, proves that A = A, in other 
words, A is real. 

An elementary argument then shows that for all real 
values of x, y and z, R{x, y, z) — R{x, y, z) — R(x, y, z), 
implying that the coefficients of the polynomial R must 
all be real. 

This completes the proof of the existence part of The- 
orem 1. 

Let us now prove that the decomposition ([2J) is unique. 

Assume we have two decompositions 

P(x,y,z) = A LiL 2 ---Ln + Q R 

= A' L[L' 2 ---L' e + Q R'. (14) 

Our goal is to show that each linear factor L\, in the 
second decomposition occurs as a factor Li in the first 
decomposition as well, modulo a possible rescaling. A 
given line L' if — intersects the quadratic Q = in a 
pair of conjugate points p and p. Because p and p satisfy 
both L'^ = and Q = 0, they satisfy P = as well. 
Turning our attention to the first decomposition, because 
p and p satisfy both P — and Q = 0, they satisfy 
L\L 2 ■ ■ ■ Lg — as well. Hence p must satisfy one of the 
lines Li = 0, and because the line's coefficients are real, 
p must satisfy that same line. But we saw earlier that 
a pair of conjugate points p and p determines a unique 
line modulo normalization (recall the essentially unique 
solution to Equations ©)• Therefore Li is a constant 
multiple of L\, , and if the coefficients of each have been 
normalized to length one, then Li = iL^, . This proves 
the uniqueness of the factorization. 

If we evaluate the two decompositions <|14[) on the com- 
plex curve Q — we get 

A L\L 2 ■ ■ ■ Lg, = A LiL' 2 ■ ■ ■ L'g (15) 

proving that if the coefficients of the Li and the L' it are 
consistently normalized, then A — A'. It then follows 
easily that R = R' as well. 

This completes the proof that the decomposition ( |2J) 
is unique, thus completing the proof of Theorem 1. 

III. COMPUTATIONAL CONSIDERATIONS 

The proof presented in Section |n] is almost construc- 
tive, but not quite. It relies on Bezout's Theorem for the 
existence of the root pairs {pi , pT, . . . ,pg, p~[\ but does not 
say how to find them. This section fills the gap. 

The key observation is that the quadratic curve Q = 
x 2 + y 2 + z 2 = is topologically a 2-sphere. More to the 
point, it is a copy of the complex projective line CP 1 , 
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which happens to be homeomorphic to the 2-sphere. Let 
us parameterize the curve Q = as 



(i(u 2 — v 2 ), —2iuv, u 2 + v 2 ) 



(16) 



where u and v are homogeneous coordinates in CP 1 . 
Clearly the mapping (|16[1 takes all points (it, v) £ CP 1 to 
the curve Q = 0, by construction. The question is, which 
of those points happen to satisfy the given polynomial P 
as well? Write 

P(x, y, z) = P{i (u 2 - v 2 ), -2i uv, u 2 + v 2 ) (17) 

to express P as a function on CP 1 . 

If v 7^ 0, then (u,v) and (— , 1) represent the same 
point in CP 1 . If we define a = — then expression l|17|l 
effectively reduces to a polynomial in a single variable, 



P(x, y, z) = P(i {a 2 - 1), -2i a, a 2 + 1) 



(18) 



The roots of this polynomial are the desired root pairs 
{px,Px, ■ ■ ■ ,pt,pi}. 

If, on the other hand, v = 0, then (u, v) = (u, 0) ~ 
(1,0). Thus (u,v) = (1,0) may represent an additional 
root, which would not be found as a root of P(a) in (fT%|l . 

Once we have found the parameters (u, v) for all 21 
roots of P, the easiest way to group them into conju- 
gate pairs is to observe that the parameterization 1|16|) 
maps "antipodal points" (u, v), (— V, u) £ CP 1 to con- 
jugate points (x, y, z), (x, y, z) £ CP 2 . In other words, 
(a, 1) and (— l,a) <~ (— 1/a, 1) map to a pair of conju- 
gate points in CP 2 . 



IV. EXAMPLES 

To illustrate how the algorithm works in practice, let 
us apply Theorem 1 to several concrete examples. 



A. Toy Quadrupole 

Consider the quadratic polynomial 

P(x, y, z) = xy + yz + zx — x 2 — z 2 . (19) 

First dismiss the special case (u, v) = (1,0) by noting 
that the parameterization i|16[l maps (u,v) = (1,0) to 
(x, y, z) = (i, 0, 1) for which HHJ gives P(i, 0, 1) = i ^ 0. 

Now consider the general case, for which Equation l|18|l 
becomes 

ia 4 + 2(1 - i)a 3 - Aa 2 - 2(1 + i)a - i = (20) 
with roots 

a x = 1 + V2 
a 2 = 1-V2 
a 3 = i(l + V2) 

a 4 = i(l-V2). (21) 



corresponding, respectively, to the four points of CP 2 , 

Pi = (1,-1, -iV§) 

Px = (1,-1, W2) 

p 2 = (-iV2,l,-l) 

P2 = (W2,l,-1). (22) 

Solving the line equations © converts the preceding two 
pairs of conjugate points to the two lines 



Lx 
L 2 



1 x + \/kv 



= 



iy+ Jlz = 0, 



(23) 



which give us the two multipole vectors ( y i, y L, 0) and 

To find the correct A, evaluate Equation (JSJ) for, say, 
q = (l,i,0), giving 



A 



P(q) 



-l + i 



Lx(q)L 2 (q) -i + i 



(24) 



Of course any other choice for q would have given the 
same answer A = 2, just so we make sure q lies on 
the curve Q(q) — x 2 + y 2 + z 2 = and exclude q £ 

{PX,P1,P2,P1}- 

We may now write down the polynomial F from Equa- 
tion |51 namely 

F = P-XLxL 2 

— (xy + yz + zx — x 2 — z 2 ) 

= -x 2 -y 2 - z 2 , (25) 

and divide by Q = x 2 + y 2 + z 2 to get the remainder term 
R = F/Q = — 1. Thus the final decomposition promised 
by Theorem 1 becomes 

2 2 

xy + yz + zx — x — z 



(26) 



+ (x 2 +y 2 + z 2 )(-l). 

B. Toy Octopole 

The cubic polynomial 

P(x,y,z) = x 2 y + y 3 



(27) 



illustrates some non-generic behavior which may arise, 
namely the possibilities of (a) a "missing root" and (b) 
multiple roots. We will follow the same algorithm as in 
Section UV Al pointing out only the differences. 

The first difference is that the special case (u, v) = 
(1,0), corresponding to [x,y, z) — (i,0, 1), is indeed a 
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root of P in 127|l . So we record that root and proceed 
onward in search of the other roots. 

The next difference we encounter is that the polyno- 
mial 



2ia° + MoT + 2ia = 0. 



(28) 



has degree only 5, not degree 21 — 2-3 = 6 as one expects 
in the generic case. Happily, this polynomial's five roots 
supplement the one exceptional root (i,0, 1) we found in 
the previous paragraph, to give the required total of six 
roots. In other words, the unexpectedly low degree of 
the polynomial is intimately tied to the existence of the 
exceptional root (i,0, 1). 

The roots of (|28l) turn out to be {— i, i, — i, i, 0}. Un- 
like more generic polynomials, this one has multiple 
roots, implying a repeated factor in the product L1L2L3. 
Specifically, those five roots correspond to 



Pi = 


1,0) 


Pi = 


H,i,o) 


P2 = 


(+M,o) 


P2 = 


H,i,o) 


P3 = 


H,o,i), 



(29) 



and then the one exceptional root («,0, 1) completes the 
pattern 

Pi = (+»,0,1). (30) 

From here the algorithm is routine. The lines are 

L x = z = 
L 2 = z = 

L 3 = y = 0, (31) 

the scalar multiple is A = — 1, and the final factorization 
is 

x 2 y + 1/ = -l{ V ){z)(z) + (x 2 + y 2 + z 2 )(y). (32) 



quadrupole, 



Y ~* = 



Y, u = 



Vo = — i 



Y 2 



15 
32^ 
15 
87 
"5" 
Ifhr 
"5" 

167T 

~5~ 

167T 

15 

15~ 

32^ 



sm 



15 

327 {X 



sm cos ( 



(x - iy)z 



(3cos 2 6» - 1) 



(3z 2 



16tt 

z 2 )) 



(3z 2 - 1) 



/ 2 2 
(-x - y 



2z l 



sin 9 cos 9 e i<p = 



2i<p 



15 
327 



s?r (x + iy)z 
(x + iy) 2 (33) 



Using the coefficients 0,2, m for the DQ-corrected Tegmark 
map of the first-year WMAP quadrupole gives 



P(x,y,z) 



-0.01262739 x 2 + 0.02302019 xy 
0.00677625 y 2 + 0.00950698 xz 
0.01064014 yz + 0.00585113 z 2 . (34) 



Following the same algorithm, as illustrated in Sections 
II V Al and II V HI we get the polynomial 

(0.01847852 -i 0.00950698) 
-(0.04604038 + i 0.02128027) a 
-0.04065752 a 2 
+ (0.04604038 - i 0.02128027) a 3 
+ (0.01847852 + i 0.00950698) a 4 

leading to the lines 

Li = 0.939660 x + 0.187066 y + 0.286437 z 
L 2 = -0.437088 x + 0.792820 y + 0.424724 z 



(35) 



0. (36) 



Converting the coefficients of these lines to spherical co- 
ordinates gives multipole vectors 



>(2,i) _ 



>(2,2) _ 



( 11.26°, 16.64°) 
= (118.87°, 25.13°) 



(37) 



in full agreement with those that CHS found using their 
tensor algorithm (Equation 3 of Q). 

An analogous computation for the octopole yields mul- 
tipole vectors v^- 3 ' 1 ', {i( 3 ' 2 ) and £>( 3 ' 3 ', again in full agree- 
ment with those reported in Equation 3 of 3] . 



C. WMAP Quadrupole and Octopole 



Our first task here is to convert a given set of coeffi- 
cients ai m to a homogeneous harmonic polynomial. Con- 
verting the standard spherical harmonics Y" 1 to poly- 
nomials m x, y and z is easy. For example, for the 



V. HOW WELL DO THE WMAP 
QUADRUPOLE AND OCTOPOLE ALIGN? 

Following p|, we define the quadrupole plane normal 
vector 



w 2,l,2 = -(2,1) x -(2,2) 



(38) 
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and the three octopole plane normal vectors 



^3,1,2 = ^(3,1) 


X 


(,(3,2) 




w 3 ' 2 ' 3 = 


X 


0(3,3) 




^3,3,1 = -(3,3) 


X 




(39) 


Still following 3], we judj 




the alignment 


of the 



quadrupole plane with the three octopole planes via the 
dot products 



as well as completely independent computer code. This 
effectively rules out the possibility of error in computing 
the l-in-1000 estimate, forcing us to take that estimate 
quite seriously. 

We should point out that our Monte Carlo simulations 
chose random quadrupoles and octopoles independently 
of each other, relative to spherically symmetric distribu- 
tions on the spaces of all spherical harmonics of degree 2 
and 3, respectively. In other words, we used independent 
Gaussian coefficients a£ m . 



Ax = 


\w 2 ^ 2 


, w 3,2,3| 










, w 3,3,l| 






A 3 = 


\w 2 ^ 2 




(40) 


VI. AN OPEN QUESTION 



Finally, in contrast to Q (which in its preprint form con- 
tained some statistical flaws), we let the sum 



S = A 1 + A 2 



A* 



(41) 



provide a numerical measure of how well the quadrupole 
plane aligns with the octopole planes. For the DQ- 
corrected Tegmark map, the sum evaluates to Sq — 2.395. 

To judge how unusually large So is, wc evaluated S for 
100000 random quadrupoles and octopoles, and found 
that only 118 trials produced S > S . This 99.9% con- 
fidence level, while weaker than the incorrectly stated 
results of 1, is completely consistent with Huterer and 
Starkman's revised statistical analysis Q. 

Like Schwarz, Starkman, Huterer and Copi, we find 
this result astonishing. In particular we find it difficult 
to believe that the quadrupole and octopole align so well 
merely by chance. Whether the alignment is imposed by 
the global topology of a small finite universe, or whether 
it is due to some previously unknown solar system effect, 
or whether it is merely the result of an error in the col- 
lection and/or processing of the WMAP data, remains 
to be seen. 

In the meantime, we emphasize that our simulations 
use an entirely different algorithm from that of 0, 0| 



Corollary 2 motivates a broader question, first raised 
by Copi, Huterer and Starkman [8j. One would like to 
decompose an arbitrary real- valued function / : S 2 — > R, 
for example the temperature function on the microwave 

sky, as a sum / = YleLo T\i=i ^,») • ^ n °tfier words, 
this approach would bypass the spherical harmonics en- 
tirely, and instead write the function / directly as the 
sum of totally factored polynomials A 
each degree £. 

Corollary 2 almost makes such a factorization possi- 
ble. For example, if we approximate the microwave sky 
temperature by the sum of its first 837 multipoles, T — 



IL=i L t,i> one for 



E837 
£=0 L~im=- 



• ae^Y™, then Corollary 2 lets us re-express 



it as 



T = Efco nLi ■ Tne question is, what 
happens when we add in the 838 t?l spherical harmonic? 
For sure we will add an 838" 1 term As38 IIi=i -^838, i to our 
factored-polynomial decomposition. And almost surely 
the 836" 1 term will change significantly as it inherits the 
remainder R from the newly added 838" 1 term. But what 
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